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Abstract. We present a new set of stellar yields obtained from rotating stellar models at solar metallicity covering 
the massive star range (12-60 M ^)- The stellar mod els were calculated with the latest version of the Geneva 
stellar evolution code described in lHirschi et al.l i2004T) . Evolution and nucleosynthesis are in general followed up 
to silicon burning. The yields of our non-rotating models are consistent with other calculations and differences 
can be understood in the light of the treatment of convection and the rate used for 12 C(a,7) ie O. This verifies 
the accuracy of our calculations and gives a safe basis for studying the effects of rotation on the yields. The 
contributions from stellar winds and supernova explosions to the stellar yields are presented separately. We then 
add the two contributions to compute the total stellar yields. Below ~ 30 A/q , rotation increases the total metal 
yields, Z, and in particular the yields of carbon and oxygen by a factor of 1.5-2.5. As a rule of thumb, the yields 
of a rotating 20 Mq star are similar to the yields of a non-rotating 30 Mq star, at least for the light elements 
considered in this work. For very massive stars (~ 60 Mq), rotation increases the yield of helium but does not 
significantly affect the yields of heavy elements. 
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1. Introduction 

Stellar yields are a crucial input for galactic chemical evo- 
lution. It is therefore important to update them when- 
ever significant changes appear in stellar evolution mod- 
els. Recent yield calculations at solar metallicity have 
been conducted by a few groups IR ausche r et alJ 2002; 
iLimongi fc Chieffil 120031: iThielemann et all 11996^ Over 
. the last ten years, the development of the Geneva evolu- 
tion code has allowed the study of the evolution of rotat- 
ing stars until carbon burning. The models reproduce very 
well many observational fe atures at various rnetallic ities, 
like surface enrichments l|Mevnet fc Maed erl |2002 |), ra- 
tios b etween red and blue supergiants ijMaeder fc Mevnetl 
l200l[) and the population of Wo lf-R ayet (WR here i nafter ) 
stars l|Mevnet fc Maederl Eoolf) . In Ifflrschi et~aTl l|2004|h 
we described the recent modifications made to the Geneva 
code and the evolution of our rotating models until sili- 
con burning. In this paper, the goal is to calculate stellar 
yields for a large initial mass range (12-60 Mq) for rotat- 
ing stars. In Sect. 2, we briefly present the model physical 
ingredients and the calculations. In Sect. 3, we describe 
the method and the formulae used to derive the yields. 
In Sect. 4, we discuss the wind contribution to the yields. 
Then, in Sect. 5, we present our supernova (SN) yields of 
light elements calculated at the pre-supernova stage. In 
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Sect. 6, we describe and analyse the total stellar yields 
(wind + SN) and compare our results with those found in 
the literature. 

2. Description of the stellar models 

The computer model us ed to calculate the s tellar models 
is described in detail in lHirschi et alJ ((2004) . Convective 
stability is determined by the Schwarzschild criterion. 
Convection is treated as a diffusive process from oxy- 
gen burning onwards. The overshooting parameter is 
0.1 Hp for H and He-burning cores and otherwise. 
On top of the meridional circulation and secular shear, 
an additional instability induced by rotation, dynamical 
shear, is introduced in t he model. The react ion rates are 
taken from the NACRE l|Angulo et alJ ll9991 compilation 
for the experimental rates and from the NACRE web- 
site (http://pntpm.ulb.ac.be/nacre.htm) for the theoreti- 
cal ones. 

Since mass loss rates are a key ingredient for the yields 
of massive stars, we recall here the prescriptions used. The 
changes of the mass loss rates . M, with rotation are take n 
into account as explained in iMaeder fc Mevnetl l|2000ajh 
A s reference mass loss rat es, we adopt the mass loss rates 
of lVink et al.l l)2000ll200l|) who account for the occurrence 
of bi-stability limits which change the wind properties and 
mass loss rates. For the domain not cover ed by these au- 
thors we use the empirical law devised bv lde Jager et all 
(.1989). Note that this empirical law, which presents a dis- 
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Table 1. Initial abundance (in mass fraction) of the chem- 
ical elements. 



Element 


Mass fraction 


Element 


Mass fraction 


L H 


0.705 


24 Mg 


5.861D-4 


3 He 


2.915D-5 


25 Mg 


7.70D-5 


4 He 


0.275 


26 Mg 


8.84D-5 


i2 C 


3.4245D-3 


28 Si 


6.5301D-4 


13 C 


4.12D-5 


32 g 


3.9581D-4 


14 N 


1.0589D-3 


36 Ar 


7.7402D-5 




4.1D-6 


40 Ca 


5.9898D-5 


!6 


9.6195D-3 


44 Ti 





17 o 


3.9D-6 


48 Cr 





18 Q 


2.21D-5 


52 Fe 





20 Ne 


1.8222D-3 


56 Ni 





22 Ne 


1.466D-4 







continuity in the mass flux near the Humphreys-Davidson 
limit, implicitly accounts for the mass loss rates of LBV 
stars. For the non-rotating models, since the empirical 
values for the mass loss rates are based on stars cov- 
ering the whole range of rotational velocities, we must 
apply a reduction factor to the empirical rates to make 
them correspond to the non-r otating case. The same re - 
duction factor was used as in iMaeder fe Mevnetl l|200lj) . 
Dur ing the Wolf-Rayet ph ase we use the mass loss rates 
bv lNugis fc Lamer!] pOOO - ). These mass loss rates, which 
account for the clumping effects in the winds, are smaller 
by a factor of 2-3 than the mass loss rates used in our 



grids dMevnet et a 


]|l994| 


Maeder fc MevnetJ 


2000a) 



). Wind anisotropy (described in 
was only taken into account for 



M ^ 40 Mq since its effects are only important for very 
massive stars. 

The initial composition of our models is given in 
Table ^ F° r a given metallicity Z (in mass fraction), 
the initial helium mass fraction Y is given by the rela- 
tion Y — Y p + AY/AZ ■ Z, where Y p is the primordial 
helium abundance and AY/AZ the slope of the helium- 
to -metal enrichment law. W e used the same values as 
in IMaeder fc Mevnetl l)200l|) i.e. Y p = 0.23 and AY/AZ 
= 2.25. For the solar metallicity, Z = 0.02, we thus 
have X — 0.705 and Y = 0.275. For the mixture of the 
heavy elements, we adopted the same mixture as the one 
used to compute the opacity tables for solar composition. 
For elements heavier than Mg, we used the values from 
lAnders fc Grevessd l|l989j) . 

We calculated stellar models with initial masses of 12, 
15, 20, 25, 40 and 60 Mq at solar metallicity, with initial 
rotation velocities of and 300 kms -1 . The value of the 
initial velocity corresponds to an average velocity of about 
220 kms" 1 on the Main Sequence (MS) which is v ery close 
to th e average observed value (see for instance iFukudal 
Il982|) . The calculations start at the ZAMS. Except for 
the 12 Mq models, the rotating models were computed 
until the end of core silicon (Si) burning and their non- 
rotating counterparts until the end of shell Si-burning. For 



the non-rotating 12 Mq star, neon (Ne) burning starts at 
a fraction of a solar mass away from the centre but does 
not reach the centre and the calculations stop there. For 
the rotating 12 Mq star, the model ends after oxygen 
(O) burning The e volution of the models is described in 
iHirschi et all ll2004h . 



3. Yield calculations 

In this paper, we calculated separately the yield contribu- 
tions from stellar winds and the SN explosion. The wind 
contribution from a star of initial mass, to, to the stellar 
yield of an element i is: 



TOn wind 

ln Pim 



t (m) 



M(to, t)[Xf(m, t) - Xf] dt 



(1) 



where r(m) is the final age of the star, M(m, t) the mass 
loss rate when the age of the star is equal to t, Xf(m,t) 
the surface abundance in mass fraction of element i and 
X® its initial mass fraction (see Tabled . Mass loss occurs 
mainly during hydrogen (H) and helium (He) burnings. 
Indeed, the advanced stages of the hydrostatic evolution 
are so short in time that only a negligible amount of mass 
is lost during these phases. 

In order to calculate the SN explosion contribution to 
stellar yields of all the chemical elements, one needs to 
model the complete evolution of the star from the ZAMS 
up to and including the SN explosion. However, elements 
lighter than n eon are marginally modified by explosive nu- 
cleosy nthesis ijChieffi fc Limongill2003t iThielemann et all 
1996) and are mainly determined by the hydrostatic evo- 
lution while elements between neon and silicon are pro- 
duced both hydrostatically and explosively. In this work, 
we calculate SN yields at the end of core Si-burning. We 
therefore present these yields as pre-SN yields. The pre- 
SN contribution from a star of initial mass, m, to the 
stellar yield of an element i is: 



pre-SN 
m Pirn 



m(r) 



m(rcm) 



[Xi{m r ) -Xf]dm r 



(2) 



where m(rem) is the remnant mass, m(r) the final stellar 
mass, X® the initial abundance in mass fraction of element 
i and Xi(rn r ) the final abundance in mass fraction at the 
lagrangian mass coordinate, m r . 

The remnant mass in Eq. |3 corresponds to the fi- 
nal baryonic remnant mass that includes fallback that 
may occur after the SN explosion. The exact determina- 
tion of the remnant mass would again require the sim- 
ulation of the core collapse and SN explosion, which 
is not within the scope of this paper. Even if we had 
done the simulation, the remnant mass would still be a 
free parameter because most explosion models still strug- 
gle to reproduce explosions l|janka et al.ll2003l and refer- 
ences therein). Neverth eless, the latest multi-D simula- 
tions l|janka et alJl2004|) show that low modes of the con- 
vective instability may help produce an explosion. When 
comparing models to observations, the remnant mass is 
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Table 2. Initial mass (column 1) and initial rotation ve- 
locity [kms -1 ] (2), final mass (3), masses of the helium 
(4), carbon-oxygen (5) cores, the remnant mass (6) and 
lifetimes [Myr] (7) for solar metallicity models. All masses 
are in solar mass units. An "A" in the second column 
means that wind anisotropy was taken into account. 



Mini 


Vi a i 


Mflnal 


M a 


Mco 




Tlifc 


12 





11.524 


3.141 


1.803 


1.342 


18.01 


12 


300 


10.199 


3.877 


2.258 


1.462 


21.89 


15 





13.232 


4.211 


2.441 


1.510 


12.84 


15 


300 


10.316 


5.677 


3.756 


1.849 


15.64 


20 





15.694 


6.265 


4.134 


1.945 


8.93 


20 


300 


8.763 


8.654 


6.590 


2.566 


10.96 


25 





16.002 


8.498 


6.272 


2.486 


7.32 


25 


300 


10.042 


10.042 


8.630 


3.058 


8.67 


40 





13.967 


13.967 


12.699 


4.021 


5.05 


40 


300A 


12.646 


12.646 


11.989 


3.853 


5.97 


60 





14.524 


14.524 


13.891 


4.303 


4.02 


60 


300A 


14.574 


14.574 


13.955 


4.323 


4.69 



usually chosen so that the amount of radioactive 56 Ni 
ejected by the star corresponds to the value determined 
from the observed light curve. So far, mostly ID models 
are used for explosive nucleosynth esis but a few groups 
have developed multi-D models fsee lTravaglio et alT 2004: 
iMaeda fc Nomotoll2003fl . Multi-D effects like mixing and 
asymmetry might play a role in determining the mass cut 
if some iron group elements are mixed with oxygen- or 
silicon-rich layers. 

In this work, we used the relation betwe en Mqq and 
the remnant mass described in lMaederHl992ft . The result- 
ing remnant mass as well as the major characteristics of 
the models are given in Tabled] The determin ation of M a 
and Mco is described in lHirschi et alJ l)2004|) . 

We do not follow 22 Ne after He-burning and have to 
apply a special criterion to calculate its pre-SN yield. 
During He-burning, 22 Ne is produced by 18 0(a,7) and 
destroyed by an a-captures which create 25 Mg or 26 Mg. 
22 Ne is totally destroyed by C-burning. We therefore con- 
sider 22 Ne abundance to be null inside of the C-burning 
shell. Numerically, this is done when the mass fraction of 
4 He is less than 10~ 4 and that of 12 C is less than 0.1. 

Once both the wind and pre-SN contributions are cal- 
culated, the total stellar yield of an element i from a star 
of initial mass, m, is: 



mp, 



pre-SN 



™?C nd 



(3) 



m Pim corresponds to the amount of element i newly syn- 
thesised and ejected by a star of initial mass m (see 
lMaederlll99l . 

Other authors give their results in ejected masses, EM: 



EM,, 



(m) rm(r) 

MXi(m,t)dt+ Xi(m r )dm r (4) 

J m(rem) 



and production factors (PF) fsee lWooslev et ai1 ll995): 

f im = EM im /[Xf{m - m(rem))] (5) 



We also give our final results as ejected masses in order to 
compare our results with the recent literature. 

4. Stellar wind contribution 

Before we discuss the wind contribution to the stellar 
yields, it is instructive to look at the final masses given 
in Table H (see also Fig. 16 in iHirschi et~ai]l2004[) . We 
see that, below 25 Mq, the rotating models lose signif- 
icantly more mass. This is due to the fact that rotation 
enhances mass loss and favours the evolution into the red 
supergiant phase at an early stage during the core He- 
burning phase (see for example iMaeder fc Mevnetl2000b|) . 
For WR stars (M > 30Mr^), the new prescription by 
iNugis fc Lamer si l|2000|) . including the effects of clumping 
in the winds, results in mass loss rates that ar e a fac- 
tor of two to three smaller than the rates from lLangerl 
l|l989|) . The final masses of very massive stars (> 60 Mq) 
are therefore never small enough to produce neutron stars. 
We therefore expect the same outcome (BH formation) for 
the very massive stars as for the stars with masses around 
40 Mq at solar metallicity. 

The wind contribution to the stellar yields is presented 
in Table The H-burning products (main elements are 
4 He and 14 N) are ejected by stellar winds in the entire 
massive star range. Nevertheless, in absolute value, the 
quantities ejected by very massive stars (M > 40 Mq) are 
much larger. These yields are larger in rotating models. 
This is due to both the increase of mixing and mass loss 
by rotation. For M < 40 Mq, the dominant effect is the 
diffusion of H-burning products in the envelope of the star 
due to rotational mixing. For more massive stars (M > 
40 Mq), the mass loss effect is dominant. 

The He-burning products are produced deeper in the 
star. They are therefore ejected only by WR star winds. 
Since the new mass loss rates are reduced by a factor of 
two to three (see Sect. EJ, the yields from the winds in 
12 C are much smaller for the prese nt WR s t ellar models 
compared to the results obtained in iMaederl <ll992l) . As is 
shown below, the pre-SN contribution to the yields of 12 C 
are larger in the present calculation and, as a matter of 
fact, t he new 12 C total yields are larger than in Maederl 
l|l992|) . In general, the yields for rotating stars are larger 
than for non-rotating ones due to the extra mass loss and 
mixing. For very massive stars (M > 60 Mq), the situ- 
ation is reversed for He-burning products because of the 
different mass loss history. Indeed, rotating stars enter into 
the WR regime in the course of the main sequence (MS). 
In particular, the long time spent in th e WNL phase (WN 
star showing hydrogen at its surface, iMevnet fc Maederl 
2003) results in the ejection of large amounts of H-burning 
products. Therefore, the star enters the WC phase with 
a smaller total mass and fewer He-burning products are 
ejected by winds (the mass loss being proportional to the 
actual mass of the star) . 

Since 16 is produced even deeper in the star, the 
present contribution by winds to this yield are even 
smaller. 12 C constituting the largest fraction of ejected 
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Table 3. Initial mass and velocity and stellar wind contribution to the yields (mp™ ld ). All masses and yields 
are in solar mass units and velocities are in kms -1 . "A" in column 1 means wind anisotropy has been included in the 
model. Z is the total metal content and is defined by: Z— 1 — Ai H — A.3 Hc — A4 He . 
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1.54E-4 


1.02E-3 


-2.87E-4 


9.20E-7 


-2 92E-6 


-3.54E-7 


() 


o 


15, 


300 


-5.23E-5 


3.25E-1 


-6.73E-3 


5.78E-4 


1.39E-2 


-7.63E-3 


1.63E-5 


-3.64E-5 


-9.37E-7 








20, 





-1.06E-5 


5.21E-2 


-1.21E-3 


1.93E-4 


2.46E-3 


-1.43E-3 


1.18E-6 


-6.28E-6 


-8.61E-7 








20, 


300 


-1.56E-4 


1.27E+0 


-1.73E-2 


1.22E-3 


4.30E-2 


-2.75E-2 


2.03E-5 


-1.01E-4 


-2.25E-6 








25, 





-6.02E-5 


3.95E-1 


-6.39E-3 


4.38E-4 


1.64E-2 


-1.07E-2 


3.39E-6 


-3.88E-5 


-1.80E-6 








25, 


300 


-2.47E-4 


2.97E+0 


-2.52E-2 


1.22E-3 


7.94E-2 


-5.48E-2 


1.03E-5 


-1.68E-4 


-2.99E-6 


2.72E-4 





40, 





-4.16E-4 


4.65E+0 


-4.48E-2 


7.11E-4 


1.45E-1 


-1.07E-1 


-1.10E-5 


-3.04E-4 


-5.18E-6 








40, 


300A 


-5.83E-4 


7.97E+0 


1.60E+0 


5.42E-4 


1.73E-1 


3.34E-1 


-3.93E-5 


-4.11E-4 


1.76E-6 


7.35E-2 


2.18 


60, 





-8.33E-4 


9.41E+0 


2.52E+0 


4.15E-4 


2.37E-1 


4.45E-1 


-6.24E-5 


-6.32E-4 


2.46E-6 


1.28E-1 


3.33 


60, 


300A 


-1.07E-3 


1.52E+1 


3.01E+0 


3.12E-4 


3.09E-1 


3.99E-1 


-9.85E-5 


-8.10E-4 


3.39E-6 


1.67E-1 


3.89 



Table 5. Pre SN contribution to the yields 

( m Pim _SN ) of solar metallicity models. Continuation of 
Table 0] Note that 20 Ne yields are an upper limit and 
may be reduced by Ne-explosive burning and that 24 Mg 
yields may also be modified by neon and oxygen explosive 
burnings. See discussion in Sect. 16.31 



M ini 


^ini 


( 20 Ne) 


22 Ne 


( 24 Mg) 


12 





1.05E-1 


6.80E-3 


6.75E-3 


12 


300 


1.58E-1 


1.55E-2 


1.36E-2 


15 





1.10E-1 


1.60E-2 


6.06E-2 


15 


300 


2.26E-1 


3.33E-2 


4.27E-2 


20 





4.83E-1 


3.64E-2 


1.28E-1 


20 


300 


6.80E-1 


4.26E-2 


1.19E-1 


25 





8.41E-1 


5.22E-2 


1.38E-1 


25 


300 


1.08E4-0 


2.24E-2 


1.48E-1 


40 





1.36E+0 


5.58E-3 


1.52E-1 


40 


300A 


1.42E4-0 


6.93E-3 


1.20E-1 


60 





1.81E+0 


6.58E-3 


1.73E-1 


60 


300A 


1.75E+0 


7.60E-3 


1.44E-1 



metals, the conclusion for the wind contribution to the 
total metallic yield, Z, is the same as for 12 C. 

5. Pre-SN contribution 

As said above, our pre-SN yields, mp?™~ N , were calcu- 
lated at the end of core Si-burning using the remnant 
mass, A/rcm, given in Table [3 We therefore concentrate 
on yields of light elements which depend mainly on the 
evolution prior to core Si-burning. 

Before discussing the pre-SN yields, it is interesting to 
look at the abundance profiles at the pre-SN stage pre- 
sented in Fig. and at the size of helium and carbon- 
oxygen cores given in Tabled The core sizes are clearly 
increased due to rotational mixing. We also see that as 
the initial mass of the model increases, the core masses 
get closer and closer to the final mass of the star. M a 
reaches the final mass of the star when M- m i > 40 Mq for 



non-rotating models and when Mini ^ 25 Mq for rotat- 
ing models. Mco becomes close to the final mass for both 
rotating and non-rotating models for Mi„; > 40Af Q . 

The pre-SN yields are presented in Tables0]andEl One 
surprising result in Tabled is the negative pre-SN yields 
of 4 He (and of 14 N) for WR stars. This is simply due to 
the definition of stellar yields, in which the initial compo- 
sition is deducted from the final one. As said above, Mco 
becomes close to the final mass for M; ni > 40 Mq. Since 
the CO core is free of helium, it is then understandable 
that the pre-SN yields of 4 He for WR stars is negative. 

5.1. Carbon, oxygen and metallic yields 

If mixing is dominant (M < 30 Mq), the larger the initial 
mass, the larger the metallic yields (because the various 
cores become larger). Rotation increases the core sizes by 
extra mixing and therefore the total metallic yields are 
larger for rotating models. Overshooting also plays a role 
in the core sizes. The larger the overshooting parameter, 
the larger the cores and the larger the yields. If we compare 
our rotating and non-rotating models, we see that the 
pre-SN total metallic yields and 12 C and 16 O yields in 
particular are enhanced by rotation by a factor 1.5-2.5 
below 30 M . 

For very massive stars (M > 30AfQ), the higher the 
mass loss, the smaller the final mass and the total metallic 
yields. The same explanations work well in general for 
carbon and oxygen. 

6. Total stellar yields 

6.1. Comparison between the wind and pre-SN 
contributions 

The total stellar yields, mp\% = mp^ SN + mp™ d (to 
be used for c hemical evolution models using Eq. 2 from 
lMaederlll992T) . are presented in Tables El and □ What is 
the relative importance of the wind and pre-SN contribu- 
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M = 60 M„, v- -=0 kmr 1 , Z = 0.02, a„, =0.1 
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m r [M ] 



6 B 10 

m r [M ] 



Fig. 1. Abundance profile at the end of core silicon burning for the non-rotating (left) and rotating (right) 15 (top), 
25 (middle) and 60 (bottom) Mq models. 
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Table 4. Initial mass and velocity and pre SN contribution to the yields (mp\^~ ) of solar metallicity models. 
All masses and yields are in solar mass units and velocities are in kms" 1 . "A" in column 2 means wind anisotropy 
has been included in the model. Z is the total metal content and is defined by: Z= 1 — Ai H — A3 Hc — A4 Ho . 
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Table 6. Initial mass and velocity and total stellar yields (mp^~ N + mp™ d ) of solar metallicity models. All 
masses and yields are in solar mass units and velocities are in kms -1 . "A" means wind anisotropy has been included 
in the model. Z is the total metal content and is defi ned by: Z= 1 — Ai H — A3 Hc — Xi Hc . These are the yields to be 
used for chemical evolution models using Eq. 2 from l|Maederlll992]) . 
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tions? Figure displays the total stellar yields divided by 
the initial mass of the star, as a function of its initial 
mass, m, for the non-rotating (left) and rotating (right) 
models. The different shaded areas correspond from top to 
bottom to for 4 He, 12 C, 16 and the rest of the heavy 
elements. The fraction of the star locked in the remnant 
as well as the expected explosion type are shown at the 
bottom. The dotted areas show the wind contribution for 
4 He, 12 C and 16 0. 

For 4 He (and other H-burning products like 14 N), the 
wind contribution increases with mass and dominates for 
M > 22Mq for rotating stars and M > 35A/ Q for non- 
rotating stars, i. e. for the stars which enter the WR stage. 
As said earlier, for very massive stars, the SN contribution 
is negative and this is why p^ em is smaller than p7^ n - 
In order to eject He-burning products, a star must not 
only become a WR star but must also become a WC star. 
Therefore for 12 C, the wind contributions only start to be 
significant above the following approximative mass lim- 
its: 30 and 45 Mq for rotating and non-rotating mod- 



els respectively. Above these mass limits, the contribu- 
tion from the wind and the pre-SN are of similar impor- 
tance. Since at solar met allicity, no WO star is produced 
l|Mevnet fc Maederl2004|) . for 16 0, as for heavier elements, 
the wind contribution remains very small. 

6.2. Comparison between rotating and non-rotating 
models 

For H-burning products, the yields of the rotating models 
are usually higher than those of non-rotating models. This 
is due to larger cores and larger mass loss. Nevertheless, 
between about 15 and 25 Mq, the rotating yields are 
smaller. This is due to the fact that the winds do not expel 
many H-burning products yet and more of these products 
are burnt later in the pre-supernova evolution (giving neg- 
ative SN yields) . Above 40 Mq , rotation clearly increases 
the yields of 4 He. 

Concerning He-burning products, below 30 Mq , most 
of the 12 C comes for the pre-SN contribution. In this mass 
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Fig. 2. Stellar yields divided by the initial mass, p\^, as a function of the initial mass for the non-rotating (left) and 
rotating (right) models at solar metallicity. The different total yields (divided by m) are shown as piled up on top 
of each other and are not overlapping. 4 He yields are delimited by the top solid and long dashed lines (top shaded 
area), 12 C yields by the long dashed and short-long dashed lines, 16 yields by the short-long dashed and dotted- 
dashed lines and the rest of metals by the dotted-dashed and bottom solid lines. The bottom solid line also represents 
the mass of the remnant (M™^/ m) . The corresponding SN explosion type is also given. The wind contributions are 
superimposed on these total yields for the same elements between their bottom limit and the dotted line above it. 
Dotted areas help quantify the fraction of the total yields due to winds. Note that for 4 He, the total yields is smaller 
than the wind yields due to n egative SN yie lds (see text). Preliminary results for masses equal to 9, 85 and 120 M Q 
were used in this diagram (see 



:gative blN yiel 
lHirschil E004). 



range, rotating models having larger cores also have larger 
yields (factor 1.5-2.5). We notice a similar dependence on 
the initial mass for the yields of non-rotating models as for 
the yields of rotating models, but shifted to higher masses. 
Above 30 Mq, where mass loss dominates, the yields from 
the rotating models are closer to those of the non-rotating 
models. The situation for 16 and metallic yields is similar 
to carbon. Therefore 12 C, le O and the total metallic 
yields are larger for our rotating models compared 
to our non rotating ones by a factor 1.5—2.5 below 
30 M Q . 

Figure 01 presents the stellar yields convolved with 
the Salpeter initial mass function (IMF) (dN/dM oc 
M -2 ' 35 ). This reduces the importance of the very massive 
stars. Nevertheless, the differences between rotating and 
non-rotating models remain significant, especially around 

20 Mq. 

6.3. Comparison with the literature 
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solar metallicity models. Continuation of Table H3 Note 
that 20 Ne yields are an upper limit and may be reduced 
by Ne-explosive burning and that 24 Mg yields can also 
be modified by neon and oxygen explosive burnings. See 
discussion in Sect. 
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We compare here the yields of the non-rotating mod- 
els with other authors. For this purpose, the ejected 

masses, EM, defined by Eq. g| in Sect. 02 are presented iThielemann et alJ(ll996l) CTNH96l lRauscher et"alll|2002t 

in Tables El and El Figure [j sh ows the compari son with fRHW02) and IWooslev fc Weaver! l|l995|) fWW95). For 

four other calculations: iLimongi fc Chieffil i|2003r (LC03), LC03, we chose the remnant masses that are closest to 
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Fig. 4. Left: total ejected masses (EM) of 1 H, 4 He, 12 C and 16 as a function of the initial mass for different non- 
rotating models at solar metallicity. Right: total ejected masses (EM) of 16 0, 20 Ne and 24 Mg as a function of the CO 
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ours (models 15D, 20B, 25A). The uncertainties related Therefore, before we compare our results with other mod- 
to convection and the 12 C(a, 7) 16 reaction are dominant. 
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Table 8. Initial mass and velocity, remnant mass and total ejected masses (EM) of solar metallicity models. All 
masses are in solar mass units and velocities are in kms -1 . "A" in column 2 means wind anisotropy has been included 
in the model. Note that this table is given for comparison with other recent publications and does not correspond to 
our definition of yields (see Sect. [31 - 
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els, we briefly mention here which treatment of convection 
and 12 C(a,7) 16 rate other authors use: 



— iLimongi k. Chieffil l)2003|) use Schwarzschild criterion 
(except for the H convective shell that forms at the 
end of core H-burning, where Ledoux criterion is used) 
for convection witho ut overshooting. Fo r 12 C(a, 7) 16 0, 
they use the rate of iKunz et alJ l|2002() (K02). 

— iThielemann et aP l)l996|) use Schwarzschild criterion 
for convection witho ut overshooting. For 12 C(a, 7) 16 0, 
they use the rate oflCaughlan et all i|l985l) (CF85). 

— IWooslev fc Weaver! l)l995l) use Ledoux criterion for 
convection with semiconvection. They use a relatively 
large diffusion coefficient to model semiconvection. 
Moreover non-convective zones immediately adjacent 
to convective regions are slowly mixed over the or- 
der of a radiation diffusion time scale to approxi- 
mately allow for the effects of conv ective overshoot. For 
12 C(q , 7) 16 0, they use the rate of lCaughlan fc Fowled 
dl988l) fCF88) multip lied by 1.7. 

— iRauscher et all ( 2002^ use Ledoux criterion for convec- 
tion with semiconvection. They use the same method 
as WW95 for s emiconvection. For 12 C(a,7) 16 0, they 
use the rate of iBuchmannl l|l996j) (BU96) multiplied 
by 1.2. 

— In this paper (HMM04), we use Schwarzschild criterion 
for convection with overshooting. Fo r 12 C(a, 7) 16 0, we 
use the rate of lAngulo eiTa l. (1999) (NACRE). 

A comparison of the different re action rates and treat - 
ment of convection is presented in Irlirschi et all i2004\ . 
The comparison of the ejected masses is shown in Fig. 
Hfor masses between 15 and 25 M Q . The 4 He and 16 
yields are larger when respectively the helium and carbon- 
oxygen cores are larger. This can be seen by comparing 
our models with those of RHW02 and LC03 (Fig. H and 
respective tables of core masses). 

For 12 C yields, the situation is more complex because 
the larger the cores, the larger the central temperature 



and the more efficient the 12 C(a,7) 16 reaction. If we 
only consider the effect of this reaction we have that the 
larger the rate, the smaller the 12 C abundance at the end 
of Hc-burning and the smaller the corresponding yield 
(and the larger the 16 yield). This can be seen in Fig. 0] 
by comparing our 12 C and 16 yields with those of LC03 
(we both use Schwarzschild criterion). Indeed the NACRE 
rate is larger than the K02 one so our 12 C yield is smaller. 
THN96 (who also use Schwar zschild criterion) using the 
rate of ICaughlan et all l)l985|) which is even larger, ob- 
tain an even smaller 12 C yield. When both the convec- 
tion treatment and the 12 C(a,7) 16 rate are different, 
the comparison becomes more complicated. Nevertheless, 
within the model uncertainties, the yields of various mod- 
els agree. In fact, the uncertainties are reduced when we 
use the CO core mass instead of the initial mass in order 
to compare the results of different groups. Fig. 0] (right) 
shows the small uncertainty for 16 O in relation to the CO 
core mass. This confirms the relation Mco~yields( 16 0) 
and shows that this relation holds for models of different 
groups and for models of non-rotating and rotating stars. 

We calculated the pre-SN yields at the end of Si- 
burning. Therefore, the yields of 20 Ne and 24 Mg may still 
be affected by explosive neon and oxygen burnings. 20 Ne 
yields are upper limits due to the possible destruction of 
this element by explosive Ne-burning. Figure 01 (right) 
shows that our results lies above the results of other groups 
but that the difference is as small as differences between 
the results of the other groups. 24 Mg yields are also close to 
the results of other groups who included explosive burn- 
ings in their calculations. By comparing our results for 
20 Ne and 24 Mg with the other groups mentioned above, we 
see that the difference between our results and the results 
of other groups is as small as the differences betw een the 
19, 20 and 21 M Q models o flRauscher et al.lll2002i) an d dif- 
ferences between for exa mple IRauscher et alJ (|2002^ and 
ILimongi fc Chieffil l)2f)f)3h . This means that our yields for 
these two elements are good approximations even though 
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Table 9. Total ejected masses (EM) of solar metallic- 
ity models. Continuation of Table [S] 
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explosive burning was not followed in this calculation. For 
24 Mg, it is interesting to note that rotation increases sig- 
nificantly the yields only for the 12 Mq models and that, 
in general, rotation slightly decreases the 24 Mg yields in 
the massive star range (see Table[7|and Fig. fright). This 
point is interesting for chemical evolution of galaxies since 
it goes in the same d irection as observational constraints 
l|Francois et al.ll2004^ . 

For 17 yields, all recent calculations agree rather 
well and differ from the WW95 results because of the 
change in the react ion rates (especially 17 0(p, a) 14 N, see 
lAubert etai1ll996|) . 18 and 22 Ne are produced by en- 
captures on 14 N. As said in Sect. 22 Ne is not followed 
during the advanced stages and we had to use a special 
calculation for its yield. Our 22 Ne v alues are neve rtheless 
very close to other calculations (see lHirschll2004j) . 

7. Conclusion 

We calculated a new set of stellar yields of rotating stars 
at solar metallicity covering the massive star range (12-60 
Mq). We used for this purpose the latest version of the 
Genev a stellar evolution code described in iHirschi et alJ 
l)2004() . We present the separate contribution to stellar 
yields by winds and supernova explosion. For the wind 
contribution, our rotating models have larger yields than 
the non-rotating ones because of the extra mass loss and 
mixing due to rotation. For the SN yields, we followed 
the evolution and nucleosynthesis until core silicon burn- 
ing. Since we did not model the SN explosion and explo- 
sive nucleosynthesis, we present pre-SN yields for elements 
lighter than 28 Si, which depend mostly on the evolution 
prior to Si-burning. Our results for the non-rotating mod- 
els correspond very well to other calculations and differ- 
ences can be understood in the light of the treatment of 
convection and the rate used for 12 C(a, 7) 16 0. This as- 
sesses the accuracy of our calculations and assures a safe 
basis for the yields of our rotating models. 

For the pre-SN yields and for masses below ~ 30 Mq, 
rotating models have larger yields. The 12 C and 16 O yields 



are increased by a factor of 1.5-2.5 by rotation in the 
present calculation. When we add the two contributions, 
the yields of most heavy elements are larger for rotat- 
ing models below ~ 30 Mq. Rotation increases the total 
metallic yields by a factor of 1.5-2.5. As a rule of thumb, 
the yields of a rotating 20 M© star are similar to the yields 
of a non-rotating 30 Mq star, at least for the light ele- 
ments considered in this work. When mass loss is dom- 
inant (above ~ 30A/q) our rotating and non-rotating 
models give similar yields for heavy elements. Only the 
yields of H-burning products are increased by rotation in 
the very massive star range. 
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